Scattering in abrupt heterostructures using a position dependent mass Hamiltonian 
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Transmission probabilities of the scattering problem with a position dependent mass are studied. 
After sketching the basis of the theory, within the context of the Schrodinger equation for spatially 
varying effective mass, the simplest problem, namely, tranmission through a square well potential 
with a position dependent mass barrier is studied and its novel properties are obtained. The solutions 
presented here may be adventageous in the design of semiconductor devices. 
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INTRODUCTION 



One dimensional quantum wells (QW) and their analysis have played an increasingly significant role in various 
applications as well as the understanding of the properties of a variety of semiconductor devicesQ, 0, 0, 0- The 
motivation for studying these problems is the recent developments in the nanofabrication of semiconductor devices, 
where one observes QW with very thin layersQ- The effective mass of an electron(hole) in the thin layered QW varies 
with the composition rate. In such systems, the mass of the electron may change with the composition rate which 
depends on the position. Therefore, the corresponding Schrodinger equation should be formulated in a correct form. 

Exact and quasi-exact solvability of the position dependent mass (PDM) Schrodinger equation has been the subject 
of recent interest El B , IE IE El El E^ El E)< It provides a useful model for the description of many physical 
systems [H El lljj llSl [19| . Although it has been solved for a number of potentials and masses, the general solution 
has not yet been completed for square well potentials. Here we suggest a model that has been easily related to the 
QW structures with various PDM models. We will demonstrate a number of promising applications of the model. 

Potential device applications, as well as purely scientific interest, provide the motivation for studies of the nature 
of the transport properties of the PDM electron through the barriers or wells. For realistic transport properties in 
semiconductors, the usual Schrodinger equation has to be replaced by the more general equation [24]: 

Q {m a pmPpm? + rn 1 pm 13 pm a ) + V(z) - E^j ip(z) = (1) 

with the constraint over the parameters: a + f3 + j = — 1. In applications, the spatial variation of m is either neglected, 
or, alternatively various special cases of Q have been suggested in the literature HHHI^. ^ n this article we focus 
on abrupt heterostructures. It has been proven 

El 

that for sharp heterostructures a — 7; otherwise the wavefunction 
is forced to vanish at the heterojunction boundary which is clearly an unphysical result. 

In contrast to the solution of the PDM Schrodinger equation including Coulomb, Morse, harmonic oscillator, etc. 
type potentials the study of the PDM Schrodinger equation including a constant potential has not 

attracted much attention in the literature. Such quantum systems have been found to be useful in the study of 
electronic properties of semiconductors. Generally, analysis of the scattering problem with PDM is based on the 
investigation of the simple problems, and it was pointed out that the transmission probability no longer tends to 
unity when incoming energy goes to infinity. The fundamental question remains open: whether the behavior of the 
transmission probability is generic or if it depends on the properties of the mass. To answer this question one has to 
obtain a general expression for the transmission probability by solving l(T)l. for an arbitrary mass. 

The rest of the paper is organized as follows. In section 2, we outline a specific formulation of the exactly solvable 
PDM Schrodinger equation to derive a general expression for the transmission amplitude of the wave through the 
square barrier. In section 3, we apply our model to calculate the transmission coefficient of the wave through the 
barrier for various spatially varying effective masses. Finally, a summary of the work and conclusions are drawn in 
section 4. 
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THEORY 



A typical QW structure is composed of a semiconductor thin film embedded between two semi-infinite semiconductor 
materials. For a compositional QW, the well material can be generated by alternate deposition of thin layers. For 
example, in a GaAs/A£ x Gai- x As QW there exists a wide GaAs well, followed by an Al x Ga\- x As barrier and a 
GaAs narrow well. The mole fraction x varies along the z-axis, therefore the mass of the electron may vary along the 
z-axis. The simplest model of the QW is that of a step potential and mass, both showing discontinuities at the same 
given point and constant inside and outside the well. Here we suggest a model by taking into account the spatial 
variation of the mass inside the barrier or well. Let us consider a potential barrier of width, d. The structure may be 
generated by continuously changing the alloy composition x of Al x Ga\- x As from x — to x = 0.32. The relation 
between alloy composition x and coordinate z is given bv[28l 

0.32z 2 

X = — • < 2 > 

Now we turn our attention to the PDM Schrodinger equation As we mentioned before, the continuity condition 
forces a = 7 = and /?=—!. With these choices the PDM Schrodinger equation QJ, takes the form: 



P—P + V a -Ej tjj{z) = 0, d > z > (3) 

where Vq is the constant potential associated with the barrier height, and E is the energy of the particle. In spite of 
its simple appearance the Schrodinger equation Q cannot be solved analytically for arbitrary to. We note here that 
an exact solution of including a constant potential can be obtained when a — 7 = — 1/4 and [3 — —1/2, but in 
this case continuity conditions can not be satisfied. We look instead at the problem from a different point of view. 
Instead of the potential Vq let us introduce the following potential |22T ]. 

r(,) = % + ^(m»-*£), d>z>0 (4) 

where to is a function of z and ml and m" denote first and second derivatives of to with respect to z. At this point it 
is worth mentioning that we will be interested in the potential which has a less pronounced cusp. Now, the potential 
resembles a square barrier or well with smooth walls. The additional term is small compared with the original potential 
Vo and does not change the shape of the potential. It is obvious that the conditions are satisfied for smoothly varying 
mass. With the potential (0} the Schrodinger equation can be exactly solved with a simple coordinate transformation 
and the wave function is given by 

iP(z) = (Cie- lfc/(z) + C 2 e lfc/W ) mi (5) 

where the function f(z) is defined as f(z) — J yfmdz and k — ^(E — Vq). 

The results given above can easily be used to solve the Schrodinger equation including well and/or barrier potentials. 
Let us illustrate our procedure on a simple example. Consider the potential barrier 

> z, z > d 

V{z) = { T/ _ , h * f m „ ? m >*\ d>2>0 ( 6 ) 




0" - 7 -£ 



with mass barrier 



[too 0>z, z>d 
m[ ' \ m{z) d> z > 1 > 

We assume that the mass of the particle too is constant outside the barrier. Mass of the particle inside the barrier 
to(z) is an arbitrary function of z. The general solution of the Schrodinger equation yields: 

r A ie lk ' z + A 2 e- ik ' z z < 

iP(z) = } {A ? ,e- lkf ^ + A A e ikf ^) mi d > z > (8) 
I A b e lk ' z z> d 
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where k' — E / h 2 and A4 are constants. For an abrupt heterostructure the continuity conditions are given bv|l8| 



m a ip(z) — continuous, vn? ^—m a i\){z) — continuous. (9) 

dz 



The transmission coefficient, T, and reflection coefficient, R, are defined by 

T=^! T + R=l (10) 

l^il 2 l^il 2 

Using elementary quantum mechanical methods, algebraic computation applying boundary conditions, will lead to 
the following expression which is related with the transmission coefficient: 

M = 

64fcfc'm m(0) 7 /4 m (d)5/4j/( d ) 

64fcfc'm m(0) 7 /4 TO (d)5/4/'((f) 
and the coefficient related with the reflection of the wave: 

~A~ Y ~ K* + {d)K*_{d)e 2lk ^ d ) - K*_{{))KX{d)e 2lk fW 



(11) 



(12) 



where K± are given by 



K±(a) = [4fc'm(a) 2 ± 4km m(a) f (a) - im m'(a)] (13) 



K±(a) is conjugate of K±(a). The transmission and reflection coefficients can be computed using the relations © 
through 111 21 . In the following section we will illustrate our model using some explicit examples. 

EXAMPLES 

In this section we discuss the dependence of the transmission probability on the position dependent mass by various 
choices of the mass m(z). We give several examples for systems with different position dependent masses. Our 
criterion for the selection of masses is that the shape of the original potential does not change and the square root of 
m(z) is analytically integrable. Moreover we made an attempt to include mass functions that are frequently used in 
the literature. In order to demonstrate our procedure, let us begin by considering the following spatially dependent 
effective masses found to be useful for studying transport properties in semiconductors: 

m a {z) = mo(<r + 5z 2 ) 

m b (z) = moae^' (14) 
m c {z) — mo(<7 + tanh(\/^z)) 
fd + 5z 2 



md{z) = m 



1 + Sz 2 



where 6 is the length scale parameter and a is a dimensionless parameter. Through out this section the parameters 
are chosen a = 0.0665, 6 = 0.0835, and Vb = lOOme^, height of the barrier and width of the barrier d — 100 A. 

It can be seen from figure the potential (0J closely resembles a square barrier with smooth walls for the masses 
m a {z), mb{z) and m c (z). We remark that when the mass rapidly changes with position z, the shape of the potential 
profile has a pronounced cusp. The potential J2J which includes the rapidly changing mass function m,d(z) can be 
plotted as shown in figure ^ We explicitly calculate transmission probability of the scattering problem employing 
various physically meaningful spatially varying effective masses in the following. 
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FIG. 1: Effect of the position dependent mass on the potential profile. The long dashed line, dashed line and dotted line show 
the effect of the m a ,m c , and raj, respectively, on the potential profile. The change in the potential profile due to rrib(z) is 
plotted with long dashed lines and it is negligible. 
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FIG. 2: The transmission coefficient T for a potential barrier with height Vb and a mass discontinuities. 



Mass barrier: m(z) = ami 

Consider now a simple mass barrier such that the mass changes at the potential discontinuities, but inside and 
outside the barrier it is a constant. In this case the potential Vb remains the same. Since the tunnelling effect is not 
qualitatively modified by the mass discontinuity, we have to leave aside the case where E < Vq. In the case E > Vq 
the calculation for transmission coefficient can easily be done from the relation and a plot of the transmission 
probability is illustrated in figure[5]for various mass ratios. In the plot we defined the quantities: 

TOi ifh E 

=a > d = i TT ' UJ = TT ( 15 ) 

ra V m oVo u 

The graph shows clearly for too > Tni the transmission coefficient no longer tends to unity when E goes to infinity, 
but it becomes an oscillating function of E, as is discussed in^J |2lJ]. In figure the curve denoted by a = 1, 
corresponds to the plot of transmission coefficient for m = 0.0665too. This is the conduction band edge effective mass 
of the electron in the Go. As structure. In the following we compute the transmission coefficients for the mass functions 
given in (|14H . 
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FIG. 3: The transmission coefficient T for a potential barrier with height Vo and various position dependent mass discontinuities. 



Mass Barriers: m a (z) , mt(z) , m c (z) and mj(z) 



The mass functions given in l|14ll are used in various fields of physics. We mention here that mass function m a (z), 
may be useful to analyze the structures GaAs/ Ai x Ga\- x As. For example the effective band mass of the electron in 
the barrier can be writtenjisl l3ll] as 



m(x) = mo (0.0665 + 0.0835a;) 



(16) 



The relation between alloy composition x and the coordinate z is given in For comparison we calculated trans- 
mission coefficients by using the relations HlOf) through (|12(1 and they are illustrated in figure 



CONCLUSION 



In summary, we have discussed the exact solvability of the PDM Schrodinger equation including a constant potential. 
We have recovered a general expression for the transmission coefficient of the wave through the square potential barrier. 
We have presented calculations of transmission coefficients for various spatially varying effective masses. 

Within the framework of the effective mass approximation, in some previous works 32] the electron was assumed 
to be confined in a square infinitely high potential well. In fact, a finite height potential well model is more realistic 
for describing the motion of the electron in the QW[33l 0| . It is obvious that the model described in this article can 
easily be modified to study QW structures and super lattices |35) . 
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